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Abstract: In the Hammersley-Aldous-Diaconis process infinitely many particles sit in 
M and at most one particle is allowed at each position. A particle at x, whose nearest 
^ . neighbor to the right is at y, jumps at rate y — x to a position uniformly distributed in the 
^ ! interval {x,y). The basic coupling between trajectories with different initial configuration 
induces a process with different classes of particles. We show that the invariant measures 
' for the two-class process can be obtained as follows. First, a stationary M/M/1 queue 
P_( . is constructed as a function of two homogeneous Poisson processes, the arrivals with rate 
! A and the (attempted) services with rate p > X. Then put the first class particles at 
the instants of departures (effective services) and second class particles at the instants of 
unused services. The procedure is generalized for the n-class case by using n — 1 queues 
in tandem with n — 1 priority-types of customers. A multi-line process is introduced; it 
consists of a coupling (different from Liggett's basic coupling), having as invariant measure 
J> ' the product of Poisson processes. The definition of the multi-line process involves the dual 
points of the space-time Poisson process used in the graphical construction of the system. 
^ . The coupled process is a transformation of the multi-line process and its invariant measure 
'vh I the transformation described above of the product measure. 

^ ■ Keywords and phrases. Multi-class Hammersley-Aldous-Diaconis process, multiclass 

! queuing system, invariant measures 
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1. The Hammersley-Aldous-Diaconis process 

The state space X is an appropriate subset of locally finite subsets of M. Elements of X 
are called configurations and elements of a configuration rj E X are called particles. The 
Hammersley-Aldous-Diaconis (^had ) process (Hammersley [15], Aldous and Diaconis [1]) 
can be informally described by saying that a particle sitting at r G M waits an exponentially 
distributed random time with rate equal to the distance to the nearest particle to its right, 
located at r' > r (say) to jump to a site uniformly distributed in the interval [r, r']. 
Alternatively, bells ring at space-time points at rate 1 and when a bell rings at (t", t), the 
nearest particle to the left of r at time t— jumps to r. 

A Harris graphical construction of the process is the following. Let G f2 be a homoge- 
neous rate-1 Poisson process in the space-time space M x (or later in M^), where Vt is the 
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set of locally finite subsets of M^. We shall use points to refer to space-time elements of u. 
If 7] is the particle configuration at time t— and (r, t) G u, then at time t the configuration 
jumps to ?7 \ {u} U {r}, where u = u{rj, r) is the nearest particle in rj to the left of r. This 
construction is well defined in a finite region as the points can be well ordered by time [1] . 
The construction in M was performed by Aldous and Diaconis [1] and then by Seppalainen 
[2T] in the state space X = {q : \\ms~>oo [Oj / ^ = oo} that we adopt. Here | ■ | counts 
the number of elements of a finite set. Homogeneous Poisson processes in M give mass 1 
to X. 

For fixed initial configuration rj and points cu, the process (?7t, t > 0) is a deterministic 
function of rj and uj denoted rjt = $(t,?7,a;), t >0. In this case we say that the process is 
governed by u with initial configuration rj. It satisfies 

r]t = ^(t - s,r]s,rsuj), (1) 

for all < s < t, where r is the time translation operator defined by TgU = {{x,t — 
s), e u}. 
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Figure 1. Harris construction. Stars represent the space-time Poisson 
points in u. Balls represent particle positions at times and t. 



Coupled and multiclass process. A joint construction of HAD processes {{rjl, t > 0), i = 
l,...,n) with initial configurations r/'^,...,?]" governed by the same points u is called 
coupled process. It is defined by 

Vl = ^t,v\u;) (2) 

We slightly abuse notation writing t] = (77^, . . . , 77") and rjt = <I>(t, 77, ci;) G A"". In the 
coupled process when a point is at the closest particle to the left of x in each 

marginal jumps to x simultaneously. 

Consider initial ordered particle configurations rj^ C ■ ■ ■ C 77" and run the coupled 
process. The order is maintained at later times. Define = {{rj^, . . . ,7/") ^ X : C 
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■ ■ ■ C r/"}. Define R : X"^^ ^ Af" by 

{Rvf = v' \ v'-' (3) 

The process C,t = Rvt is called the multiclass process. Particles in Q are called i-class 
particles. The multiclass process is just a convention to describe a coupled process with 
ordered initial configurations. The map R is invertible; its inverse is given by {R~^^)^ — 
^1 U ■ ■ ■ U The process governed by uj with initial configuration ^ is defined by 

it = i:{t,i,u) = R^t,R-^i,u). (4) 

Invariant measures. The Poisson process with density A is an invariant measure for the 
HAD process for all A > [T]. Our main result is to construct invariant measures for the 
multiclass process. The resulting measure coincides with the law of the departure process 
of a stationary multiclass-customer queue system. 

Let A and S be particle configurations in X . Think of M as time and construct a 
continuous time random walk (Z^,, r G M) jumping one unit up at times in A and one unit 
down at times in S. We fix Zq = but since we are interested in the increments, the 
position at a given time is not important. The increments satisfy 

Zr- Zs = \An[r,s)\-\Sn[r,s)\ (5) 

Assume that — >■ Too as r — > ±oo; this implies Zr visits each site a finite number of 
times. Let S) be the times in S that attains a new record down and D{A^ S) its 
complement: 

U{A,S) := {reS:Zr< inf Z,/} 

r'<r 

D{A,S) ■= S\U{A,S) (6) 
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Figure 2. Arrivals and services. Each arrival is linked to its departure time, 
running the queue according to the FIFO (first-in-first-out) schedule. 
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We interpret A as arrival times of customers to a one-server queue and S as tlie times 
at wliicli service occurs. At tliese service times tliere are two possibilities: eitfier tfiere is a 
customer in tlie system, producing a departure, or tliere is no customer, in this case there 
is an unused service. We collect the departure times in D and the unused service times 
in U. If A and 5" are stationary Poisson processes with rates A < A', respectively, then Zr 
satisfies the conditions above and D{A, S) are the departure times of a stationary M/M/1 
queue {Qr , r G M) defined by = {Z^ — inf^/^^ Z,.')^; here Qr is the number of customers 
in the queue at time r. This construction of {Qr, r G M) is standard (reference ??). The 
process is stationary Markov and satisfies 

Qr - Qr- = l{r eA}-l{reD : > 0} (7) 

The operator D has the following properties: (i) D{a^,a^) C a^. (ii) If C 
then D{a^,a'^) C D{a^,a'^). (iii) If and are independent one-dimensional Poisson 
processes of densities and with < p^, then D{a^,a'^) also is a one-dimensional 
Poisson process with density p^. This is just Burke's Theorem for a M/M/1 queue [3]. 

Let be the law of a Poisson process of rate A. For p^ < ■ ■ ■ < p" define u : = 
i/p X ■ ■ ■ X . Let a = {a^, . . . , a") G A"" be a multi-line configuration with law u. Define 
a sequence of operators Z^*-"^ : X as follows. Let D^^\a^) = a^, and then recursively 

for n >2, let 

= D (^("-^^(aS . . . , a""^), a") . (8) 

The configuration D^'^\a^, a^, . . . , a") represents the departure process from a system of 
(n — 1) queues in tandem. The arrival process to the first queue is . The service process 
of the /cth queue is a^^^ , for k = 1, . . . ,n — 1. Finally, for k = 2, . . . ,n — 1, the arrival 
process to the kth queue is given by the departure process of the {k — l)st queue. This is 
known system of -/M/l queues in tandem. 

Note D^'^\a^,a'^) = D{a^,a'^). By applying (i)-(iii) above repeatedly, we obtain that 
D^'^\a^, . . . , a") C Z)("~^)(q;^, . . . , a") C ■ ■ ■ C a" and if a^, . . . , a" are independent one- 
dimensional Poisson processes of densities < ■ ■ ■ < p", then D^'^\a^ , . . . , a"') also is a 
one-dimensional Poisson process with density p^. 

Define the configuration t] = [rj^, . . . , 77") by 

r/'= = D("-'=+^)(a^«'^+^...,«"). (9) 

By construction rj'' C rj^'^^ for all /c = 1, ...,?7, — 1 and for each k, rf^ has marginal 
distribution . 

Define the map C : X"" 1— > X"-^ by Ca = rj. Define the measure vr on X"-^ as the law 
of rj (that is, vr = Cu). Define the multiclass measure p as the law of ^ = Rrj, that is, 
p = RCu. Call M = RC. See Figure O 

Theorem 1. If a has law v with p^ < ■ ■ ■ < p", then the law Mu of ^ = Ma is invariant 
for the multiclass process C,t defined in Equivalently, the law Cv of Ca is invariant for 
the coupled process rjt defined in ([^. 
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Figure 3. Construction of a coupled configuration rj and a multiclass con- 
figuration ^ from a multi-line configuration a. C,'' particles are colored ac- 
cording to its class. In the middle picture, arrivals have been linked with 
their departure times in each server, according to the FIFO schedule and the 
customer class. 

To prove Theorem [T] we introduce later a new process at = {aj, . . . , a") G A"" called 
multi-line process. Proposition H] shows that the product measure u is invariant for at and 
Proposition [S] shows that if a is the initial configuration for the multi-line process at then 
{Cat, t > 0) is the coupled process with initial configuration Ca. These propositions imply 
that Cu is invariant for the multiclass process. 

The right density of a configuration rj E X is the limit as r goes to infinity of the 
number of //-particles in [0, r] divided by r, if this limit exists. Analogously define the left 
density and call simply the density of rj the result when the left and right density are the 
same. The density of an ergodic measure in A* is the unique value A such that the set of 
configurations with density A has u probability one. 

Uniqueness. The next result says that the measure Cu is the unique invariant measure 
for the coupled process in the set of ergodic measures on Af" with marginal densities p. Its 
domain of attraction includes the ergodic measures with densities p. Let p^ < ■ ■ ■ < p", 
p = (p^, . . . ,p") and p' = {p^,p^ — p^ . . . ,p" — p"~^). Let M.^{p) be the set of ergodic 
measures on A"" such that the ki\i marginal has density p^ for k = 1, . . . , n. 

Theorem 2. Assume the conditions of TheoremUi Then, 

1 ) The measure Mv is the unique invariant measure for the multiclass process in Ai"'{p') . 
The multiclass process starting with a measure in 7W"(p') converges weakly to Mv as t ^ 
oo. 

2) The measure Cv is the unique invariant measure for the coupled process in M.'^i^pY 
The coupled process starting with a measure in 7W"(p) converges weakly to Cu as t ^ oo. 
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For all A > and almost all realizations u; of the points, there is a unique stationary 
realization of the had process at density A governed by uj. This is shown by the next 
theorem; notice that the notation rjt in this theorem is used for the HAD process on X, 
that is, with n = 1. 

Theorem 3. For each A > there exists an essentially unique function H\ mapping 
elements uo of Vt to had trajectories {r]t,t G M) such that: 

(i) The induced law of {rjt.t G M) = Hxioj) is stationary in time. 

(ii) The marginal law of rjt for each t is space- ergodic with particle density A. 

(iii) With probability 1, {rit,t ^M.) is a HAD evolution governed by uo. 

(Here "essentially unique" means that if H[ is another function satisfying the three con- 
ditions, then Hx{uj) = H'^{uj) with probability 1). Then in fact the marginal law of rjt for 
each t is . 

To stress the dependence on A call rj^ the process constructed in Theorem [31 The 
construction implies rj^ G 7]^ if A < A'. The union of rj^ in A is a countable dense set of M. 
It consists on the space coordinates of the points in u with time coordinate less than t. 

Theorems [2] and [3] are proven in Section |H based on Proposition [10] which considers 
the coupled process {t]l,'>]'f) starting with two independent configurations with the same 
density A. If one of the configurations is a Poisson process and the other comes from an 
ergodic distribution, then, with probability one, the density c(t) of positions where the 
two configurations differ at time t is deterministic and converges to zero as t grows. The 
proof of the proposition follows an argument of Ekhaus and Gray [7] as implemented by 
Mountford and Prabhakar [T^ . 

Multi-line process. The stationary realization (77^, t G M) of the had process governed 
by u and density A of Theorem [3] induces a point configuration on M x M consisting on 
the space-time positions of the particles just before jumps. Gator and Groeneboom [1] 
call them dual points and prove that they have the same Poisson law as u for all A > 0. 
The proof is based on two facts: (a) the reverse had process with respect to z/^ is also 
a had process, with jumps to the left instead of to the right and (b) a trajectory of the 
process uniquely determines the points governing it. The dual points and a new density 
produce a new stationary realization of the had process, which in turn produces new dual 
points. Repeating the procedure n — 1 times, we construct a stationary realization of the 
multi-line process with n lines. Fix densities p^, . . . ,p" and a realization of the points u. 
Gall (a", t G M) the stationary realization of the had process at density p" governed by 
u! (which is well-defined by Theorem [3]), and, for A; = n — 1, . . . , 1 set u;''" = dual points 
generated by p'^ and u^^^ with = uj (it is convenient for later reasons to start at the 
bottom) and a^~^ as the stationary realization of the HAD process at density p^~^ governed 
by uj^. Define (at, t G M) as the multi-line (stationary) process given by at = {aj, . . . , a"). 
The generator of the multi-line process is given later in (fTSj) . The multi-line process can 
be constructed for times in [0, 00) for any initial distribution; we do so in Section [2l In this 



7 



ffi- 



©■ 



© A 



©- 



© * 



©-* 



time 



© ± 



© ± 



time t 



LO, a\ and u)"^ 



Figure 4. The dual points are represented by circles 



case the n — 1 first marginals do not follow the had dynamics unless the initial distribution 
is a product of Poisson processes. We prove two results in Section [2) 

Proposition 4. Let p = [p^, . . . , p^) be positive densities. The product of Poisson processes 

1 k 

X ■ ■ ■ X pP is the unique invariant measure for the multi-line process with marginal 
densities p. 

Proposition 5. Take a multi-line configuration a such that the multiclass configuration 
Ma is well defined. Let at he the multi-line process with initial configuration a governed by 
uj. Then (Mctt, t>Q) is the multiclass process with initial distribution Ma and {Cat., t > 
0) is the coupled process with initial distribution Ca, both governed by uo. 

The proof of Theorem [1] is a consequence of Proposition H] (without the uniqueness part) 
and Proposition O 

Prahofer and Spohn [20] and Patrik Ferrari [13] introduce a multilayer dynamics based 
on the dual points and related to our multi-line process by a rotation of 45° (roughly 
speaking). See the pictures in [TS] . 

Multiclass Burke's Theorem. The invariant measure for the multiclass had is also a 
fixed point for a multiclass ./M/1 queue. 

To make this more precise we need to introduce more notation in the construction of rj 
as a function C oi a. As in [11] we label multiclass configurations as a function of the 
multi-line configuration a. Let ^^'^ = a^ and define 

^^.'^ := U{a^-^,a^) (10) 

for k = 1, . . . ,n and i = 1, . . . ,k. Interpret these configurations as the service times of 
customers of k classes in the kth system; at times in ^^-^^^ customers of class i arrive to 
the kth system and are served at times in The service schedule respects the classes: 



8 



first class customers are served first, second class customers are served when there are no 
first class customers waiting, etc. End of services occur at times in a'^. The unused service 
times of the kth system are in One can think that the kth system has infinitely many 
customers of class k waiting, so that they are served at the service times unused by the 
lower class clients. With this interpretation ^^'^ = are the service times of a system with 
infinitely many customers in queue at all times. Assuming FIFO schedule (first-in-first- 
out) for customers of the same class we can link each time in a'' with the corresponding 
service time in a'^'^^. This is the meaning of the links in Figure [31 In each line k, the 
particles linked to some particle in the first line will be first class particles; particles linked 
to particles in the second line but not to particles in the first line are second class particles 
and so on. 

Let M : A"" A"" be the map between multi-line configuration a and the multiclass 
configuration ^ = (^'^'^, . . . ,^"'") using and (ITOll : 

Ma:=iC'\---,C'n (11) 

where A*" C X"^ is the subspace of X where M is well defined. This coincides with the 
map M = RC defined just after 

The nth multiclass stationary queuing system (Q", r G M) in {0, 1, . . . jg defined 

as a deterministic function of the arrivals ■ ■ ■ ,^^~^'"'~^) and the services a"^ by 

Q^U) - Q^r-U) = Hr G C-'''} - l{r e : g^(^) = for ^ < j; Q^j) > 0}{12) 

where Q"(j) is the number of customers of class j in the nth system at time r. The 
existence of a process satisfying f[T^ can be proved by induction on n by observing that at 
times r in ^"'''^ the queue is empty: Q^U) = for all j < n. The process is not Markov 
unless the arrival process is product of homogeneous Poisson and the service process is also 
Poisson. The departure times of class j in system n is given (equivalently to (fTOj) ) by 

^"J = {r G a" : g^(z) = for z < j and Q^(j) > 0} (13) 

for j = 1, . . . , n — 1. As before ^"'^ is just the set of unused service times in the nth system. 

Burke theorem [3] says that the departures of a stationary M/M/1 queue have the same 
Poisson law as the arrivals. It applies to our case for n = 2: Q^(l) is an M/M/1 queue 
with arrivals Poisson of rate and departures The extension to the multiclass 
system says that the departures of the (n — 1) classes of the nth system (.^"'^, . . . ,^"''^~^) 
have the same law as the arrivals (^"~^'^, . . . ,C"~^''^~^) to the same system. This is one of 
the consequences of the multiclass Burke's theorem as follows. 

Theorem 6 (Multiclass Burke). Fix k > 1 and for n > k let a & A"" have law v with 
< ■ ■ ■ < p"' . Then the law of {$,^'^, . . . , ^"''^') (= first k coordinates of Ma ) is independent 
ofn. 

Proof. Let 77^"^ be the coupled process with initial invariant distribution Cv^'^\ where i/^"'' 
is the product measure in with marginals v'^ , p^ < ■ ■ ■ < p^ . Since the evolution of 
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the first k coordinates is Markovian, the marginal law of the first k coordinates of 77^ 
under Cz/*-"-' is invariant for the first k coordinates of the process 7]^""^ By uniqueness of 
the invariant measure of Theorem [21 the law of these marginals must coincide with Cu^''\ 

(k) 

the invariant measure for rji , the process with k lines and the same densities. Since 
(^"■'^, . . . = first k coordinates of Rr]^^\ its law coincides with the law of Rt]^''^ = 

Ma^^\ where a^'^^ is the multi-line configuration with law z/^'^). □ 

The process in a cycle. Similar results can be proven for the multiclass process in a cycle 
Mtv The space state is the set of finite configurations contained in [0, A^]. The points uj 
are restricted to Mjv x I^^- When the configuration at time t— is r]t_ = rj and there is a point 
at [x, t), if there are no 77 particles to the left of x, then the rightmost particle of 77 jumps 
to X. The coupled process is defined as in ([2]). Take two finite particle configurations A]\f 
and Sn in Mat with \Sn\ > \Ai\f\. Extend these configurations periodically to two infinite 
configurations A and 5*. Construct {Zt, t G M), a periodic process satisfying ([5]). The 
resulting periodic configurations of departures D{A,S) and unused services U{A,S) are 
periodic and induce configurations and Un in Mat- We construct a (unique) queue 
Qr, r G Mat such that it satisfies ([7]) which has value Qr = for r e Un; actually Qr so 
constructed is the minimal process satisfying ([7]). 

Given a multi-line configuration a in XJ^ we construct a multiclass configuration ^ = 
M^a G XJ^, where Mn is defined as the map M substituting D and U by and Un- The 
analogous to Theorems [1], [2] and [6] hold for this process. Letting i = [i^, . . . , £^), the set 
A'"(p) of Theorem [2] must be substituted by XJ^{Q, the set of configurations with exactly 

particles in the kth cycle. The proof of the analogous to Theorem [2] is easy in this case. 

Regeneration properties of the invariant measure. In the two-class invariant mea- 
sure the second class particles are regeneration events. More precisely, let /i on X"^ 
be a translation invariant measure with marginal densities p^ip^- The Palm measure 
"conditioned to have a second class particle at the origin" is the measure fi defined by 

r(i^2^i fi^rC)-: fo^^ any measurable bounded set / C M, where |J| 
is the Lebesgue measure of / and 9^ is translation by r. The multiclass invariant measure 
Mv of Theorem [T] conditioned to have a second class particle at the origin satisfies: (a) 
the configuration to the left of the origin is independent of the configuration to its right; 
(b) both the distribution of first class particles to the right of the origin and the positions 
of the first plus second class particles to the left of the origin are Poisson processes (or 
product measures in the discrete-space case). These properties were proven by [5] for the 
two-class invariant measure for the totally asymmetric simple exclusion process (tasep); 
alternative probabilistic proofs can be found in [9]. Angel [2] shows that these properties 
follow easily from the two-class representation ^ = Ma that he called "collapsing proce- 
dure". In Section Owe show that as in the discrete case [TT], there are regeneration strings 
for the multiclass invariant measure Mv. 

Microscopic shocks. The process with two classes of particles is a crucial tool to define 
shock measures [inilH]- Taking a configuration (^^, ^^) distributed according to the invariant 
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measure Mu for the two-class process conditioned to have a ^"^ particle at the origin and 
considering the configuration consisting on all particles (of any class) to the left of the 
origin and the first class particles to the right of it, then the distribution of the resulting 
configuration is invariant for the process as seen from an isolated second class particle. We 
discuss this item in Section |6l 

Multiclass HAD process in Z and other discrete processes. The approach has been 
implemented for a discrete-space multiclass had process and other discrete-space processes 
on Z. In fact, this research started with the extension of Angel [2] result to the multiclass 
process for the discrete HAD and TASEP we performed in [11] and [12]. There are two main 
points in this paper: one is to show that the approach works also in continuous space (the 
HAD process is the most natural model to study tasep like questions in R) ; the other is to 
show some results that are only announced in [12]. In particular we show the uniqueness 
of the multiclass invariant measure with given marginal densities. The use of dual points 
is more intuitive in the continuous had than in the discrete processes. Indeed, the dual 
points in this case are identified by the trajectory of the process, while in the discrete cases 
it is necessary to appeal auxiliary spin-fiip processes to identify the dual points. 

The the existence of an invariant measure for the two-classes TASEP is first show by 
Liggett [TB] [IZ] and then computed by Derrida, Janowsky and Lebowitz [5], see also 
[22] and [1]. Angel [2] and Duchi and Schaeffer [B] description of the two-class invariant 
measures was the starting point of our multiclass version [11] and |12j . 



2. The multi-line process 

The multi-line process has configurations at = {aj, . . . , a") in A"". Given a multi-line 
configuration a and a position x G M let be the position of the closest particle of a" 
to the left of x and inductively for lines k = n — 1, . . . , 1 let be the position of the 
closest particle to the left of x^~^^. When the configuration at time t is a and there is 
an u point in (x,t), the a^-particle at jumps to x at time t and simultaneously the 
particle located at x'' jumps to x'^"^^ for k = n — 1, . . . ,1. The function J : A"" x M ^ Af" 
defined by 

(J(a,x))'= = a^U{x'=+^}\{x^} (14) 

maps the multi-line configuration a before the jumps produced by x to the configuration 
after the jumps. The generator of the process is given by 

£/(«)= / rfx[/(J(«,x)) -/(«)] (15) 

JR 

This definition is equivalent to the one given in the introduction. For t such that (x, t) is 
in u denote x"^^(t) = x and x'^(t) the position at time t— of the particle in the kth line 
jumping due to the point in (x, t). Then u'' = {(x'^(t),t) : (x, t) G u} are the dual points 
of cu^"'"^ which in turn govern the process a^. 
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We define another multi-line had process a^: in A*" and show that it is the reverse of at 
with respect to the product of Poisson processes z/. For the points govern the first line 
but producing jumps to the left: when a point appears at y, it calls the closest particle 
to the right of y, located at y^. Simultaneously the closest particle to the right of y^, 
located at y^ jumps to y^, and so on. Calling yo = y, the positions {yo, ■ ■ ■ , y-n) are defined 
as a function of a and y. The multi-line configuration obtained after the jumps produced 
by a point at y is called J*{a,y); its kth line is given by 

ir{a,r))' = a'u{y'~'}\{y'} (16) 
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FlGURE 5. Effect of a point at x for the multi-line process and the reverse 
jump (due to a point at y in the reverse). 



The operators J and J* are the inverse of each other in the following (equivalent) senses: 
J*( J(a, x), x^) = a; J{J*{a, y),y'') = a (17) 
(in the figure, = y and y"' = x.) The generator of is given by 

£*/(«)= / dy[f{r{a,y))-f{a)] (18) 

Proposition 7. For any choice of positive densities p = (p^, . . . , p"), the multi-line process 
with generator C* is the reverse of the multi-line process with generator C with respect to 
the product measure v on A"" with densities p. As a consequence, v is invariant for both 
C and C*. 

Proof. The proof is based on the observations: a) the configurations a and J{a,r), the 
configuration after a jump produced by a Poisson point at r, have the same probability 
weight; b) the rate to jump from a to J(a, r) is the same as the rate in the reverse process 
to jump from a' = J{a,r) to J*{a',y) = a; c) the rate of exiting any configuration a 
is the same for the direct and the reverse process. We formalize this with a generator 
computation. 
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It suffices to show that for bounded functions / and g depending on finite regions in 
v{gCf) = iy{fC*gy, that is, 



^{da) / g{a)dx[f{J{a,x)) — f{a)] = / ^{da) / f{a)dx[g{J*{a,x))—g{a)] 
Jr J Jr 

Since the term subtracting in both terms is the same, it suffices to prove that 

iy{da) / g{a)f{J{a,x))dx= / u{da) / g{a)f{J*{a,x))dx (19) 



By conditioning on the positions of the jumping particles (x^, . . . , x") the left term in f|T9|) 
can be written as 

I dx I ]^(pfce-^'("'^'-"')) dx\.. dx" / v^{da) g{a)f{r{a, x)) (20) 

Jr Jx^< - <x"<x J 

where is the measure z/ conditioned to have at each line k a particle at position x^ and 
no particle in {x^,x^^^). Change variables: call a' = J{a,x) so that, by (fTTj) a = J*{a',y) 
and y'^ = x'^'^^, A; = 0, . . . , n to get that fl20|) equals 



[dy [ l[{p'e-^'^y'-y'-'^)dy\..dy'^ [ uy{da')g{r{a',y))f{a') (21) 

Jr Jy<y^< -<y" j^^-^^ J 

where is the measure u conditioned to have at each line k a particle at position y'' and 
no particles in {y^~^, y^). Expression (1211) is just the right hand side of f|T9l) . □ 

The set of dual points of the had process starting with a Poisson process of rate A is a 
Poisson process on M x ]R+. This is proven in Theorem 3.1 in ^ under the title "Burke's 
theorem for Hammersley processes". We state and prove this fact in our context. 

Proposition 8. If r] ^ X has law and u is a rate-1 homogeneous Poisson process 
m M X M"*", then the dual points of the trajectory {rjt, t > 0) form a rate-1 homogeneous 
Poisson process in M. x M+. Furthermore for t > 0, rjx is independent of the dual points 
contained m M x [0,T]. 

Proof. By Proposition [7] (with n = 1), the reverse had process with respect to is 
just a HAD process with the drift to the left. Assume we have constructed the process 
iVtj < t < T) with initial law and Poisson points on M x [0, T]. The reverse process 
?7* = TjT-s-, < s < T, is governed by the dual points in M x [0, T]. Since the reverse 
process started in equilibrium z/'', the points governing the reverse process must be also 
Poisson in M X [0,T]. 

The dual points in M x [0, T] are the points governing the future of the reverse process 
?7j := rjT-t- which is in equilibrium. Hence it is a Poisson process independent of the initial 
configuration r]^ = rjT- □ 

In the proof of Theorem [3] we use Proposition [8] to construct a stationary trajectory 
(?7j, t G M) of the HAD process governed by points u with marginal law z/'*'. The dual points 
of this trajectory are, as u, a homogeneous Poisson process in M^. 
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3. Multi-line and multiclass processes 



In this section we prove that the projection of the multi-hne process on the multiclass 
space using the map M is just the multiclass process. The analogous statement is true 
for the coupled process. Recall that given an initial multi-line configuration a and a 
homogeneous rate-1 Poisson process in M x the process (a^, t > 0), where at = 
ip{t,a,u!) is a realization of the multi-line process. On the other hand, given an initial 
configuration ^ and the same u, the process (^t, t > 0) with = ^{t, C,, u) is the multiclass 
process. 

Proposition 9. Assume Ca is well defined. Then almost surely: 

C{^{t,a,uj)) = ^{t,Ca,uj) (22) 
M(^(t,a,cj)) = T{t,Ma,uj) (23) 

In particular, the law of Cat is the same as the law ofrjt. 



a3 ^ 5~Q ^CT^ 



o o mr~~^ Q_ 

_Q Q Q QO^ 



• • 

? • 0_Q • 

7,3 Q Q e Q_ 



before 
jump 



before 
jump 



after 
jump 



after 
jump 



Proof. Expression ([23]) follows then from the identity M = and dH). Fix x G M; 

it suffices to show 

H{Ca,x) = CJ{a,x) (24) 

where i?(?7, x) is the coupled Hammersley configuration obtained after a point at (x, t) (for 
some t) if rit- = rj and J {a, x) is the resulting multi-line configuration after a point at (x, t) 
if at- = a: its /cth coordinates satisfy 

{H{v,x))'' = r/\{^'=}U{x} 
(J(a,x))'- = a'^ \ {x''} U 

where 2;*^ is position of the closest rj'' particle to the left of x and x'' is the position of the 
particle jumping on line k in the multi-line process; here x^~^^ = x. 
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Let 7] = Ca he the coupled configuration obtained from a and denote a = J [a, x) and 
T] = H{ri,x). With this notation ( 12^ is Ca = rj. We first prove ( 124|) for n = 2 hues and 
then use this case to proceed by induction. 

The case n = 2. The second line of both processes are governed by the point at x with 
the same Hammersley rule. This implies they coincide in the second line: (Ca)^ = (?7)^. 
There is an service time at x and an q} particle at x^. This particle is served in line 2 
at time x unless there is another particle to the left of served at that time. In any 
case, an q} particle is served at x in line 2. This means {CqY has a particle at x. To see 
that {CqY has no particle at z-^ and that otherwise it is equal to rj we consider two cases 
depending whether x"^ serves or not an particle. 

1) The time at x"^ serves an particle. Hence, there is an rj^ particle at x^. Since 
there are no services between x^ and x (by definition of x^), there cannot be rj^ particles 
either and r]^ = r]^ \ {x^} U {x}. Lets check that also Ca is of this form. Let y be the 
position of the a^ particle served at x^ in line 2. Clearly y < x^. The a^ particles to the 
left of y are served in line 2 before x^ and hence its service time is not affected if x^ and x^ 
are translated to the right. This means Cq coincides with 77 at x and to the left of x. The 
a^ particles to the right of y are served after x. The displacement of the service time from 
x^ to X just changes the service time of the particle at y, leaving the other service times 
unchanged. The displacement of the a^ particle from x^ to x^ does not change the order 
or arrival of the particles neither their service times at a^, as those times are after x. This 
implies Ca = 77 in this case. 

2) The a"^ time at x^ does not serve an a^ particle. Then the a^ particle at x^ is served at 
some time z before x^. There are no a^ particles between x^ and x^. The translation of the 
a^ particle from x^ to x^ and the a^ particle from x^ to x do not change the service times 
of the a^ particles to the left of x^, which are served before x^. By the FIFO schedule, 
particles arrived after x^ are served after particles arrived before. Hence, the translations 
do not change the service times of the particles arrived after x^. 

The induction step. From the definition of C, we have 

r/f = /^("-'+'U«f,...,«r). (25) 

From the observation just before IQ, 77^ has distribution u'''^ . So we need to show that 
the RHS of (l25l) is a had trajectory governed by u. Since {a^, . . . , a^) is itself a multi- 
line process (with n — k + 1 lines) governed by uj, it is enough to show that, for any n, 
D^^^al, . . . , a^) is a had trajectory governed by u. But from the definitions of Z)^"^ and 
of the multi-line process, 

D^-\al . . . , an = D^'^ {D^^^'^al . . . , c^r^), <) 

This and the fact that {a], . . . , a"~^) is an (n — l)-line multi-line process governed by 
concludes the induction step. □ 
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4. Uniqueness of the invariant measure 



In this section we prove Theorems [2] and [3] using the following Proposition. It says 
that the marginals at time t of the coupled Hammersley process with initial independent 
marginals will coincide as t ^ oo if the first initial marginal is a Poisson process and the 
second one is an ergodic distribution with the same density as the first. 

Proposition 10. Let A > and rj^fj E be independent particle configurations with 
law , the Poisson process and i)^, an ergodic process in X with density \, respectively. 
Let {rjt^fit) he the coupled process with initial configuration {rj,'?]). Then with probability 
one, both rjt \ fjt and fjt \ rjt have a deterministic density denoted c{t) = c{t] u^) which is 
decreasing and converges to as t ^ oo. 

An immediate consequence of the proposition and the translation invariance of the law of 
the coupled process at each given time is that the probability of the event {Ar\r]t = Anfjt} 
converges to 1 as t — > cxd for any bounded interval A C M. 

Before proving the proposition we show how it implies Theorems [2] and [31 

Proof of Theorem [2]. Assume r] and fj in X^ are independent with laws fi and fl re- 
spectively and let {rit,fjt) be a 2?T,-coordinates coupled process starting with configuration 
[f], fj). Assume fi and fi are invariant for the coupled process and that the marginal law of 
both T)'' and rj'^ is a Poisson process with rate p'' for all k. Then, for bounded / depending 
on a bounded interval A: 



which tends to zero by Proposition [10] applied to each coordinate. This implies fi = fi. n 

Proof of Theorem [3]. First construct a double infinite realization and the correspond- 
ing points. Fix particles rj with law Poisson of rate A and Poisson points u;~ U uj~^, the 
subsets of points with negative and positive time coordinates, respectively. Run the pro- 
cess ?7t = $(t,?7,u;+) for t > using the points W^. Run the reverse process backwards 
V~t — ^*{^tjVj^~) starting from the same configuration using the points Here 
$*(— t, ?7, cj") = ?7, TR(u;~)), where TR(u;) = {(x, — t), {x,t) G u} are the points of 
uj reflected with respect to the line {t = 0}. Let the dual points D~{lj~ , rj) be the positions 
of the particles of 77^ just before jumps. By Proposition [HI D~{uj~, rj) is a Poisson process 
of points. Let u = D~{uj~,ri) U be the configuration consisting on the dual points of 
the reverse process for negative times and the original points for positive times. The points 
u are Poisson and govern a stationary process having configuration rj at time zero. This 
constructs simultaneously Poisson points u and a stationary trajectory governed by u. 




(26) 



(27) 



(28) 
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For uniqueness we need to show that if (f/j,) = (^s, — oo < s < oo) is a stationary 
evolution governed by u such that the time marginal fj^ is ergodic with density A then 
fis = rjs for all s. By Theorem [2] the marginal law of fjt is Poisson of parameter A for all t. 
Hence both {fjs) and (rjs) are space-time ergodic processes but the joint process {{rjs), {fjs)) 
is not necessarily ergodic. Fix a positive t and introduce an auxiliary process (r;^, s > —t) 
governed by u. The initial configuration rj'^f. is Poisson of rate A and independent of the 
configurations {ri_t,fi_t) at that time. By Proposition [TOl the density of tiqAtiq and fjoArjQ 
are both smaller than 2c{t) for all t. Hence the density of fjoArjo is bounded by 4c(t) for all 
t. Since c(t) converges to 0, this implies that r]o = f/o and by the same argument r]s = Vs 
for all s e M. □ 

Proof of Proposition IIOI We use an argument of Ekhaus and Gray [7j as developed in 
Section 2 of Mountford and Prabhakar [19]. Let = Vt^Vt, = Vt\Vt and = rjt\fjt. 
As in [19] we call them yellow, blue and red particles, respectively. For each t the process 
{iT^it^it) is ergodic and has the same density as . We want to show that the 
density of goes to zero as t oo. Label the particles of rjt and fjt as follows. Call rjtii) 
the position of the particle of rj. Initially ?7o(«) < fjQii + 1) for all i; the same for fj. 
The labels evolve in time depending on the color. At the time of an uj Poisson point at x, 
proceed as follows: 

1) If the closest left rj and fj particles are in then both of them jump to x carrying 
their labels. 

2) If the closest left rj particle is blue localized in xi < x, the closest fj particle is red 
localized in yi < xi, then the rj particle at xi and the fj particle at yi jump to x carrying 
their labels and change their color to yellow. If furthermore there are blue particles in 
yi < Xk < ■ ■ ■ < X2 < xi, the blue particle in Xi jumps to Xj+i, i = 2, . . . , k — 1 carrying 
the label and the color. 

2') The same as (2) by interchanging r] with fj and blue with red. 

3) If the closest r] particle is blue at position xi < x and the closest fj particle is yellow at 
position yi < xi, then the yellow particles in yi jump to x carrying the label and keeping 
the color. If furthermore there are blue particles in yi < Xk, ■ ■ ■ ,X2 < Xi, then the rj blue 
particle at Xi jumps to Xj+i, z = 2, . . . , A; — 1 carrying the label and the color. 

3') The same as (3) by interchanging r] with fj and blue with red. 

In other words, when a labeled particle becomes yellow it keeps the color for ever. 
Yellow particles behave as first class particles while blue and red particles do as second 
class particles. When blue and red particles coalesce they change the color to yellow; this 
happens in cases (2) and (2'). Since there is no creation of new particles, its density 
must be non increasing. If a labeled rj or fj particle overpasses another particle, then it 
becomes yellow. As a consequence blue and red particles can be overpassed but cannot 
overpass other particles. 
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Call everblue those particles in the initial configuration 77 that will be blue at all t > 
and everred those particles in fj that will be red at all times. The configuration of everblue 
particles has a translation invariant distribution but not necessarily ergodic. Let 7^^ C rj 
be the set of 77 particles that at time t will be blue. The law of 'j^ is ergodic with a 
(deterministic) density c{t), the density of c(t) decreases to a value c as t — > 00. Since 
the configuration of everblue particles is the intersection in t of 7^^, it has density c. For 
the same reason the configuration of everred fj particles has also ergodic law with density c. 

Lemma 11. Assume the conditions of Proposition UU If the density of everblue particles 
is strictly positive, then there exist an M < 00, an n < 00 and a 5 > such that for each 
t, the density of red particles W at time t satisfying 

i) there exist blue particles at time t in {W^ W + M], 

ii) there are at most n rjt particles in (W, W + M], 

iii) there are no rjt particles in {W + M, + M + 5], 
IS at least I /{AM). 

Proof. The existence of an M such that the the red particles W at time t satisfying (i) 
have at least density 1/M is proven in Lemma 3.1 of [19] for another process. This is the 
hard part of the Ekhaus-Gray argument. The proof applies here because it only uses the 
fact that red and blue particles cannot overpass, so everred and everblue particles maintain 
the order. Since rjt is a Poisson process, the density a{k) of rjt particles U having at least 
k rjt particles in the interval (f/, U + M) decreases exponentially with k. For n sufficiently 
large so that a{n) < 1/ (2M) the density of red particles at time t satisfying (i) and (ii) for 
this n is at least 1/(2M). Finally, as rjt is a Poisson process, the density h{5) of r]t particles 
U such that there is at least one rj particle in the interval {U + M, U + M + S\ is bounded 
above by 5p (the mean number). Take 5 sufficiently small such that h{5) < 1/(4M). Then 
the density of red particles W at time t satisfying (i), (ii) and such that there is another rj 
particle in {W + M, ly + M + 5] is smaller than 1/(4M) and the red particles W at time 
t satisfying (i-ii-iii) have at least density 1/(4M). □ 

Lemma 12. Under the conditions of Lemma [771 there exists a positive e' such that for all 
t the density of S,t minus the density of ^^^-^ is at least e' . 

Proof. Take integer times t and for each rjt particle U consider the following event in the 
space-time Poisson process u: 

iv) there are no u points in (t/, U + M] x [t, t + 1/2]. 

v) u has exactly n + 1 points (xi, ti), . . . , (x^+i, tn+i) in {U + M,U + M + 6] x [t, t + 1] 
and they are increasing in the time coordinate and decreasing in the space coordinate. 

This event has a positive probability e" and it is independent of the past up to t. This 
implies that the density of red particles at time t satisfying conditions (i) to (v) has some 
positive density at least e' = e"/8M. 
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The red particles at time t satisfying (i) to (v) will collide with a blue particle between 
t and t + 1. Hence the density of — is not smaller than e'. □ 

Proof of Proposition [TUl By contradiction. If we assume that the density of decreases 
to a non negative constant c, then by Lemma [12] this same density decreases by a fixed 
amount 5' > at each unit of time. □ 

5. Regeneration properties of the multiclass invariant measure 

Assume that there is a second class particle at the origin. This corresponds to an unused 
service in the queue, which in turn implies there are no customers in the queue at time 
0. Since the queue is Markovian, the future depends only on the number of customers at 
time 0. The attempted departures to the right of the origin is a Poisson process with rate 
p. We conclude that conditioned on having a second class particle at the origin, the first 
plus second class particles to the right of the origin form a Poisson process of rate p. By 
reversing time, the arrivals in the reverse process also form a Poisson process of rate A. 
The arrivals in the reverse queue are just the effective departures in the forward queue; 
that is, the first class particles. Hence, conditioned on having a second class particle at 
the origin, the first class particles to the left of the origin form a Poisson process of rate A 
and independent of the right process. 

For n > 3 classes there are not regeneration events but as in the discrete case, there are 
regeneration strings. Fix a vector of classes (cq, . . . , q) such that cq = n, q = 2 and for each 
class m G {3, . . . , n — 1} there exists a position jm < i with Cj„ = m and with Cj < m for all 
3 ^ {jm + 1, . . . , -^j. For instance c = (4, 1, 2, 3, 1, 2) qualifies for £ = 5 and js = 3. We call 
c a regeneration string. Let X^{c) be the set of configurations ^ in such that there are 
positions = Xq < ■ ■ ■ < such that G ^'^^ and fl (xj, Xj+i) = for = 1, . . . , n. In 
other words, for ^ G X^{c) there is a ^" particle at the origin and looking at the classes of 
the first i particles to the right of the origin, their classes follow the vector c. A translation 
invariant measure p on with densities conditioned on X'^{c) is the Palm 

measure p further conditioned on A'"(c); here the Palm measure is the measure conditioned 
on having a ^" particle at the origin defined by pf = (l/|/|p") / dp{^) J2re^"ni f i^^O 
any measurable bounded set / C M. As in the discrete case one can prove the following 
result 

Proposition 13. Let c be a regeneration string of length i and Xq < ■ ■ ■ < x^. Then, 
under the multiclass invariant measure Mu conditioned on X'^{c) the string ^ fl (— oo,Xo) 
is independent of {xi, 00). 

The proof is analogous to the two-classes case. One has to verify that after a regeneration 
string the queue is empty, so whatever will happen in the future (of xi) depends on the 
Poisson processes of the future, so it is independent of the past of xe. In particular the 
superposition of the classes to the right of the regeneration string U ■ ■ ■ U ^" U {xe, 00) 
under the conditioned measure of the Proposition is a Poisson process of rate p". See more 
details in 
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6. Shocks in had 



Liggett [16] introduced the coupled process to prove ergodic properties of the exclusion 
process; he called "discrepancies" what we call second class particles. The interest on 
the study of invariant measures for the coupled process was renewed when its association 
to shock measures become clear. A shock related to a process rjt is a (possibly random) 
position depending on t with the property that uniformly in time the asymptotic densities 
to the right and left of the shock are different. Wick [23] showed the existence of that 
position for the totally asymmetric zero range process with densities and A to the left 
and right of the origin, respectively. In fact this is equivalent to the position of the leftmost 
particle in the tasep with densities and A to the left and right of the origin, respectively. 
Then [10] and [8] (see also [18]) used the existence of an invariant measure for for n = 2 
to show the existence of a shock in the tasep. In a similar way Garcia [H] proved the 
existence of a shock measure for the had process. 

We show now a way to construct a shock measure for the HAD with asymptotic densities 
A > p to the left and right of the origin, respectively. Start with the invariant measure fi of 
Theorem [T] for the coupled process with two marginals, z/"^ marginal for C,^ and marginal 
for U.^^. Let fi be the measure /i conditioned to have a second class particle at the origin, 
and be a two-class configuration chosen according to fi. From construct a configuration 
r] by superposing the negative particles of C,^ and all the particles of except the one at 
the origin: 



Call S the map that transforms ^ in 77 and jj^x.p the law of rj so constructed. For the coupled 
process with initial configurations rj and rj' := rjU {0}, there is only a discrepancy for all 
times. This discrepancy behaves like a second class particle. Call Xt its position. The 
process {rit,Xt) is Markovian but the marginal process Xt is not. Let 77^ be the process 
defined as the translation by Xt of r]t. Call the operator that takes t, rj' and u in ri[. 
Then /i' is invariant for 77^. The proof is based on the fact that S commutes with the 
dynamics of ^t in the following sense: 



where $^ is the operator that transforms a time t, points u and initial configuration ^ in 
the configuration of the two-class process at time t as seen from the second class particle 
Txt^t- Here is the translation operator defined by t^A = {y — x, y E A}. The map $' 
takes a time t, points u and an initial (one-class) configuration 77 into the process as seen 
from a second class particle Txtlt- This works in the same way as for the tasep [T0| [8] 
and for the had [T3], so we omit the proof. 



7/ = [e^n(-oo,0)]ue'\{0} 



(29) 



5<D2(t,a;,0 = $'(t,^,^0 



(30) 
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